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Summary 


This paper addresses the problem of stabilization of a class of internally passive nonlinear 
dynamic systems using linear, time-invariant (LTI) passive controllers. Fundamental results 
on global asymptotic stability are obtained via Lyapunov-LaSalle method, using extensions 
of the Kalman- Yakubovich lemma. It is shown that a class of nonlinear time-invariant 
systems can be robustly stabilized by LTI controllers which are strictly positive real in the 
weak or marginal sense. 


1 Introduction 


«•(*(<)) = 


Let $R denote the set of real numbers, and 5?" the n- dimensional Euclidian space. Let 
£ 2 := £2(0, 00) denote the Lebesgue space of real- valued functions of time t which are 
square-integrable on the interval [0, 00), and let ££ denote the set of n-tuples of such 
functions. The truncation operator, Px, is defined as the function 

x(t) for t < T, 

0 for t > T. 

Let xx(t) := P T (x(t)) denote the truncation of x(f). The set £3. == £" 2 , 10 , 00) denotes the 
set of all functions x(t) whose truncations xx(t) are in ££ for all finite T. ££ e * s called the 
extension of the space ££ [1]. 

The inner product of x and y in is defined as: 

roo 

<x,y>= x T (t)y(t)dt. 

Jo 

The norm of x E ££ is defined as ||x|| = (< x,x >)2. The truncated inner product 
< x,y >t is defined as: 


f T 

< x, y >x=< x T ,yT >= / x T (t)y(t)dt. 

Jo 


The truncated norm of x(t) is defined as: ||x||t = ||xx||- 

We consider a class of nonlinear systems (denoted by £) affine in control, described by: 


x = f(x) + g(x)u 
y = h(x) + N(x)u 


( 1 ) 
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where the state vector x G 5R n , the input u and the output y G C™ e , f(x) G 8£ n , h(x) G 9? m , 
g(x) G 9£ nXm , N(x) G 5? mxm , and f,g,h,N smooth (i.e., infinitely differentiable, or C °° ) 
functions, and /( 0) = 0, h(0) = 0. 

We shall assume that the system E given by Eq. (1) is globally reachable and zero-state 
observable. These properties are defined next. 

Reachability and Observability 

Let $(f , x 0 , u) denote the state transition map for initial state x 0 at t = 0 and input u(t). A 
function cr(.) is said to belong to class K if it is strictly monotonic increasing and cr(0) = 0. 
Let B r denote the set {x G ?R n , ||x|| < r}. 

Reachability: A system E is reachable if there exist a constant r > 0 and a class K 
function a such that, for every x G B r , there exist a finite t\ > 0 and an input u, IHloo < 
cr(||x||), such that x = $(<i,0,u). If this holds Vx G 3?”, E is globally reachable. 

We shall assume in this paper that the system under consideration is globally reachable. 

Zero-State Observability - A system, E is said to be locally zero-state observable if, 
with u(t) = 0, there exists a neighborhood fl of 0 such that Vx G fl 

y = h(Q(t,x 0 , 0)) = 0 Vf> 0=>x = 0. (2) 

Also, if Q = 3?" the system is said to be zero-state observable. 

In addition, we shall introduce the following property which will be used in the stability 
proofs. 

Strong Zero-State Observability: A system E is said to be strongly zero-state observ- 
able if it is zero-state observable, and has the property that {lim^oo u(t) = 0, and lim^oo y\ t ) 
0} implies lim^oo x(f) = 0. 

2 Background 

This section contains the relevant definitions of passivity, as well as the state space char- 
acterization of passive nonlinear and linear systems. Some stability results in the existing 
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literature are also given. 


2.1 Passivity 

The concept of passivity was first introduced in the network theory literature. In the 
context of electrical networks, passivity has the implication that any single-port network 
consisting solely of resistors, capacitors, and inductors constitutes a passive system. Simi- 
lary, for mechanical systems, any spring-mass-damper system with non-negative damping 
coefficients is passive. 

Passivity can be defined in the input-output (10) sense or in the internal sense [2]. 
The 10 definition is more general and is applicable to a large class of systems including 
time- varying and infinite-dimensional systems. 

Input-Output Passivity: A system is said to be passive in the 10 sense if there exists a 
constant j3 such that 

< u, y > T +0 > 0 Vu € ££, VT > 0 (3) 


Strict IO Passivity: A system is said to be strictly passive (or input strictly passive) in 
the 10 sense if there exists a constant (3 and a constant e > 0 such that 

< u,y >t +(3 > e||u||j Vu € £™ e , VT > 0 (4) 

Internal passivity is usually defined for finite-dimensional systems, as follows. 

Internal Passivity: The system (E) is said to be internally passive if there exists a 
nonnegative storage function fj[.] such that 

<u,y> T >E[x(T)]-E[x( 0)] Vu € C™ e , VT > 0. (5) 

Strict (or input-strict) passivity in the internal sense is defined similarly. 

The difference between 10 passivity and internal passivity is that a non-negative stor- 
age function (which is a function of the system’s state vector) exists for the latter case. 
For finite-dimensional systems that are globally reachable and zero-state observable, 10 
passivity and internal passivity are equivalent. 
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We shall next present some passivity definitions for finite dimensional, linear, time- 
invariant (FDLTI) systems. 

2.2 Passive FDLTI Systems 

For finite-dimensional linear, time-invariant (FDLTI) systems, input-output passivity is 
equivalent to internal passivity of a minimal realization. The equivalence is a result of the 
Kalman- Yakubovich lemmafl]. Therefore, in the discussion of FDLTI systems, we shall 
use the term passivity to represent input/output passivity. For such systems, passivity 
is equivalent to “positive realness” of the transfer function [3], [4]. The concept of strict 
positive realness has also been defined in the literature, and is closely related to strict 
passivity. 

Let G(s) denote an m x m matrix whose elements are proper rational functions of the 
complex variable s. G(s) is said to be stable if all its elements are analytic in Re(s) > 0. 
Let the conjugate-transpose of a complex matrix H be denoted by H*. 

Definition 1: An m x m rational matrix G(s) is said to be 
positive real (PR) if 

(i) all elements of G(s) are analytic in Re(s) > 0; 

(ii) G(s) + G*(s) > 0 in Re(s) > 0; or equivalently, 

(iia) poles on the imaginary axis are simple and have nonnegative-definite residues, 

and 

(iib) G{ju) + G*(ju>) > 0 for to £ (— oo, oo). 

Given below are some definitions of strictly positive real (SPR) systems found in the 
literature. Definition 2, which represents the specialization of the general definition of strict 
passivity to stable LTI systems, is the strongest definition of strict positive realness. 

Definition 2: An m x m stable rational matrix (?(s) is said to be strictly passive if there 
exists an e > 0 such that 

G(ju>) + G*(ju>) > el for lv £ (— 00 , 00 ). 

Strictly passive systems require the system to have a relative degree of zero, which 
makes this a very restrictive class of systems. 
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Definition 3 (ref. [7]): An m x m stable rational matrix G(s ) is said to be strictly 
positive real in the weak sense (weak SPR, or WSPR) if 

G(ju>) + G*(ju) > 0 for u> € (— 00 , 00 ). 

A stronger definition of SPR was given in [5], [6], which requires certain additional 
conditions. Definitions 2 and 3 assume stable systems. The following definition, introduced 
in [8], allows the system to have poles on the imaginary axis. 

Definition 4 (ref. [8]): An m x m rational matrix G(s ) is said to be marginally strictly 
positive real (MSPR) if it is positive real, and 

G(ju) + G'(ju) > 0 for uj G (— 00 , 00 ). 

Definition 4 of MSPR differs from Definition 1 (PR) because the frequency domain 
inequality (>) has been replaced by the strict inequality (>). The difference between 
Definitions 3 and 4 is that the latter allows G(s) to have poles on the imaginary axis. Thus 
Definition 4 gives the least restrictive class of SPR systems. If G(s) is MSPR, it can be 
expressed as: G(s) = Gi(s) + £ 2 ( 5 ), where G 2 ( 5 ) is WSPR and all the poles of G\(s) (in 
the Smith-McMillan sense) are purely imaginary [8]. 

2.3 State-Space Characterization of Passive Systems 

The following theorem ([9], [2], [10]) gives the necessary and sufficient condition for a system 
to be internally passive, and is the nonlinear version of the well-known Kalman-Yakubovich 
lemma. 

Theorem 1 The nonlinear system (1) is internally passive if and only if there exists a 
non-negative function E(x) G C 1 , E( 0) = 0, and functions l(x) G 9ft*, and W(x) G 9ft* Xm 
for some integer k such that 

V r £(x)/(x) = -l T (x)l(x) 
g T {x)VE{x) = h{x) - 2W T {x)l(x) 

W T (x)W{x) = \[N(x) + N T (x)]. (6) 

4U 


5 



For FDLTI systems, Theorem 1 results in the Kalman-Yakubovich lemma. In [7], the 
Kalman- Yakubovich lemma was extended to WSPR systems, and in [8], it was extended 
to MSPR systems. These extensions are given next. ([A, £?, C, D] denotes an nth-ordei 
minimal realization of the m x m transfer function matrix G(s)). 

Lemma 1 (ref. [7]) G(s ) is WSPR if and only if there exist real matrices: P = P T > 
0, P € 5R nxn , L <E & mxn , W € such that 


a t p + pa = -l t l 

( 7 ) 

c = b t p + w t l 

( 8 ) 

w T w = d + d t 

( 9 ) 


where [ A , B, L, W] is minimal and F(s ) = W + L(sl — A) l B is minimum phase. 

Proof The proof can be found in [7]. □ 

Let [A 2 , B 2 , C 2 , D ] denote an nth-order minimal realization of G^s), the stable part ol 
G(s). The following lemma is an extension of the Kalman-Yakubovich lemma to the MSPR 
case. 

Lemma 2 (ref. [8]) If G(s) is MSPR, there exist real matrices: P = P T > 0, P (E 
9? nxn ,£ € r xnj , W G $R mxm such that Eqs. (7)-(9) hold with 

L — [0 m xni j £m X n2 ] (10) 

where [ A 2 , B 2 , £, W] is minimal and F(s) = W + L(sl — A)~ 1 B = W - f C(sl — A 2 )~ l B 2 is 
minimum-phase. 

Proof Refer to [8]. □ 

Only the sufficiency part is given in this lemma because it is of relevance in obtaining 
the stability results to be presented. Also note that in both WSPR and MSPR cases 
the storage function is given by E(x) = \x T Px, where x is the state vector of a minimal 
realization of G(s). 

2.4 Previous Stability Results 

Stability analysis of the negative feedback interconnection of two passive systems has long; 
been a topic of considerable interest. The best-known stability result for such systems 
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is the Passivity Theorem, which states that the negative feedback interconnection of a 
passive system and a strictly passive system is finite-gain (or Lf)- stable (e.g., see [1]). If 
the systems are internally passive, globally reachable, and zero-state observable, then the 
feedback interconnection is globally asymptotically stable [2]. The requirement of strict 
passivity, however, is too restrictive. In particular, a strictly passive system requires a 
direct feedthrough term in the output equation (i.e., the relative degree has to be zero). It 
is, therefore, highly desirable to weaken the requirement of strict passivity. 

To that effect, in the case of FDLTI systems, the requirement of strict passivity was 
weakened in [7] and [8]. It was shown in [8] that the negative interconnection of a PR 
system [f/(s)] and an MSPR system [(C(s)] is asymptotically stable if none of the ju;-axis 
poles of G(s) is a transmission zero of H(s). The stability holds regardless of the presence 
of unmodeled dynamics and parametric uncertainty in the systems, and is therefore robust. 
A special case occurs when G(s) is WSPR, and for that case, the resulting feedback loop 
is also asymptotically stable, as was stated in [7]. 

We shall next extend the stabilization results of [7] and [8] to nonlinear passive systems 
in a negative feedback loop with passive FDLTI systems. 

3 Stabilization of Passive Nonlinear Systems Using 
LTI Controllers 

Consider the system shown in Figure 1, which represents the negative feedback intercon- 
nection of a passive nonlinear system and an LTI system. The following theorems establish 
global asymptotic stability of the feedback system under weaker conditions, i.e., when one 
of the systems is passive, and the other is linear and WSPR or MSPR. 

Theorem 2 Suppose in the system shown in Figure 1 , E is affine in control ( see (1)), inter- 
nally passive, strongly zero-state observable, and has a radially unbounded storage function 
E(x). Then the feedback system is globally asyptotically stable if the system G(s) is WSPR. 

Proof 2 Let x denote the n-dimensional state vector of E and x denote the n-dimensional 
state vector of a minimal realization [A, B,C, D\ of G(s). Since G(s) is WSPR, from [7], 
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there exist matrices P = P T > 0, PG 3? nXn , L E 9? mXn , W E 3? mxm , such that Eqs. (7)-(9) 
axe satisfied, and [A, 5, L , W] is minimal and F(s) = W+ L(sl — A)~ l B is minimum-phase. 
Consider the candidate Lyapunov function 

V(x,x ) = E(x) + E(x) 

where E(x ) = Px is the storage function of G(s). Note that E(x) is positive definite 
because the system is zero-state observable [2]. Differentiation of V(t ) with respect to t 
and simplification using (1), and (6)-(9) leads to 

V = — d(x,u ) + y T u — z T z + y T u (11) 


where 

d{x,u) = [l(x) + W(x)u] T [l(x) + W(x)u] 

and 

z(t) = ±{Li + Wu). 

Noting that u = — y, u = y, we have 

V = —d(x,u) — z T z < 0 (12) 

i.e., V is negative semi-definite, and the closed-loop system is at least Lyapunov-stable. 
However, we will show that the closed-loop system is, in fact, globally asymptotically 
stable. 

We have established so far that V < 0. Now, V = 0 =$■ d(x(t),u(t)) = 0, and z(t) = 0. 
However, z{t) is the output (produced by the input u(t)) of the system: F(s) = -^[W + 
L(sl — A)~ l B\, which has transmission zeros only in the open left-half plane. Therefore, 
u{t) —> 0 exponentially, i.e., y(t) — > 0. Since G(s) is stable, u(t) — > 0 => y(t) — > 0. That 
is, u(t) — >■ 0 and y(t) — > 0. Since E is strongly zero-state observable and [ A , B , C, D] is 
minimal, u(t), y(t) -4 0 =$► x(t) -4 0, and x -)• 0, i.e., \\x(t)\\ -4 0, ||x(t)|| -4 0. This implies 
that V(x(t),x(t)) -4 0, i.e., it shows that V(x(t),x(t)) decreases ( V < 0) somewhere along 
the trajectories implying that V ^ 0. Therefore, such trajectories cannot exist, and by 
LaSalle’s invariance principle, the feedback system is globally asymptotically stable. □ 
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The conditions of Theorem 2 can be relaxed to allow the system G(s) to have poles 
on the imaginary axis. However, this requires the nonlinear system E to satisfy an ad- 
ditional condition. The following theorem gives the conditions under which the feedback 
interconnection of a passive system and an MSPR system is globally asymptotically stable. 


Theorem 3 Suppose in the system shown in Figure 1, E is affine in control, internally 
passive, strongly zero-state observable, and has a radially unbounded storage function E(x). 
In addition, suppose E has the property that u(t) $ £2(0,00) =£>■ limt-**, y(f) 7^ 0. Then the 
closed-loop system is globally asymptotically stable if G(s) is MSPR. 


Proof 3 Proceeding as in the proof of Theorem 2, suppose x denotes the n-dimensional 
state vector of E and x denote the n-dimensional state vector of a minimal realization 
[A, B,C, D] of G(s). Let n 2 denote the number of poles of the stable part G 2 (s) of £?($), 
and let [A 2 , B 2 , C 2 , D] denote its minimal realization. Since G(s) is MSPR, from Lemma 
2, there exist matrices P = P T > 0, P e 3R nXn ,£ € 3? mxn2 ,H r € 3R mXm , such that Eqs. 
(7)-(10) axe satisfied, and [A 2 , B 2 ,C , W] is minimal and minimum-phase. 

Consider the candidate Lyapunov function 

V{x,x) = E{: x) + E(x) 

where E(x) = ^x T Px is the storage function of G(s). Using Eqs. (6)-(10) and proceeding 
as in the proof of Theorem 2, the time derivative of V can be obtained as 

V = —d(x, u ) + y T u — z T z + y T u (13) 


where 

d(x,u ) = [l(x) + VP(x)u] T [Z(x) + W(x)u] 


and 

z(t) = -^[Li + Wu] = ±{Cx 2 + 

where x 2 is the state vector corresponding to [A 2 , B 2 ,C 2 ,D]. Noting that u = — y, u = y, 
we have 



V = — d(x,u) — z T z < 0 


(14) 
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i.e., V is negative semi-definite, and the closed-loop system is at least Lyapunov-stable 
However, we will show that the closed-loop system is, in fact, globally asymptotically 
stable. 

V = 0 =>• d(x(t),u(t)) = 0, and z(t ) = 0. However, z(t ) is the output (produced by the 
input u(t)) of the system: F(s) = ^=[ W + C(sI — A 2 ) -1 1?2], whose transmission zeros are in 
the open left-half plane. Therefore, u(t) — >■ 0 exponentially, i.e., y(t) — ► 0. Therefore, y(t) 
can consist only of exponentially decaying terms, and sinusouidal terms (including zero- 
frequency) corresponding to the jw-axis poles of G(s). If y(t) consists of any sinusoids, 
then y(t) £ 2 [0,oo), which implies that y(t) cannot go to zero; this contradicts the fact 
that y{t ) -¥ 0; therefore, y(t) can consist only of exponentially decaying terms. Since E is 
strongly zero-state observable and [A, 5,(7, D] is minimal, u(t) and y(t) — > 0 =>■ x(t ) — * 0. 
and x(t ) — ¥ 0. The rest of the proof is the same as that of Theorem 2. □ 

In Theorem 2, the requirement that u(t) £ £ 2 ( 0 , 00) =*> lim^oo y(t) ^ 0 will usually be 
satisfied if the zero dynamics of E are asymptotically stable. 

4 Concluding Remarks 

The problem of robust stabilization of a class of nonlinear passive systems was consid- 
ered. It was shown that, under the assumption of strong observability, a class of internally 
passive nonlinear systems can be stabilized by a class of linear, time-invariant passive con- 
trollers. The stability holds regardless of the presence of unmodeled modes or parametric 
uncertainties, and is therefore robust. 
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